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Interparticle interactions in charge-stabilized colloidal suspensions, of arbitrary salt concentra- 
tion, are described at the level of effective interactions in an equivalent one-component system. 
Integrating out from the partition function the degrees of freedom of all microions, and assuming 
linear response to the macroion charges, general expressions are obtained for both an effective elec- 
trostatic pair interaction and an associated microion volume energy. For macroions with hard-sphere 
cores, the effective interaction is of the DLVO screened-Coulomb form, but with a modified screening 
constant that incorporates excluded volume effects. The volume energy - a natural consequence of 
the one-component reduction - contributes to the total free energy and can significantly influence 
thermodynamic properties in the limit of low-salt concentration. As illustrations, the osmotic pres- 
sure and bulk modulus are computed and compared with recent experimental measurements for 
deionized suspensions. For macroions of sufficient charge and concentration, it is shown that the 
counterions can act to soften or destabilize colloidal crystals. 



PACS numbers: 82.70.Dd, 83.70.Hq, 05.20.Jj, 05.70.-a 



I. INTRODUCTION 

More than a century ago, it was recognized that most colloidal particles carry an electric charge (!]]. Colloidal 
macroions - typically 1 — 1000 nm in diameter - may acquire charges from surface dissociation of counterions, 
adsorption of salt ions from solution, or creation of defects in crystal lattices. Electrostatic repulsion between macroions 
suspended in a molecular fluid is one of the two chief mechanisms by which colloidal suspensions may be stabilized 
against coagulation induced by attractive van der Waals forces. 

Charge-stabilized colloidal suspensions exist in a wide variety of forms. Familiar examples include clay minerals 
(relevant to mineralogy, agriculture, and the paper industry), paints, inks, and solutions of charged micelles. Further 
examples are synthetic latex or silica microspheres which may self-assemble, if sufficiently monodisperse, into 
ordered crystals. Aside from providing valuable model systems for fundamental studies of condensed matter, colloidal 
crystals exhibit unique optical properties that have inspired a number of recent applications, e.g., nanosecond optical 
switches || , chemical sensors Q , and photonic band gap materials ^| . 

Despite the considerable and growing technological importance of charged colloids, progress in predicting macro- 
scopic properties is limited by an incomplete understanding of interparticle interactions. Most theoretical treatments of 
electrostatic interactions are rooted in the landmark theory of Derjaguin, Landau, Verwey, and Overbeek (DLVO) [||. 
The DLVO theory describes the bare Coulomb interactions between macroions as screened by the surrounding mi- 
croions (counterions and salt ions). The resulting screened-Coulomb pair potential accounts - at least qualitatively 
for a range of observed phenomena, including the dependence of coagulation rate on counterion valence and trends in 
phase stability with varying salt concentration. Recently, interest in colloidal interactions has intensified as a result 
of accumulating experimental evidence for apparent long-range attractions between macroions . 

A rigorous statistical mechanical treatment of the multi-component mixture of macroions, counterions, salt ions, 
and solvent molecules is a daunting task. Interactions in such complex systems are therefore usually treated at the 
level of effective interactions. Tracing out from the partition function statistical degrees of freedom associated with 
all but a single component, the mixture is formally mapped onto an equivalent one-component system of "pseudo- 
particles" governed by an effective state-dependent interaction |)| . Effective interactions in charge-stabilized colloids 
have been modeled by a variety of techniques, including Poisson-Boltzmann cell models [To|-[l^| density-functional 
theory Monte Carlo and Molecular Dynamics simulation p9[-p^|, and powerful ab initio methods p^ ) pq |. 

Here we adopt an alternative approach, recently proposed by Silbert et al. |^J,^5| , which exploits analogies between 
charged colloids and metals. Performing a classical trace over microion degrees of freedom and treating the electro- 
static response of the microions to the macroions within second-order perturbation theory, leads to an effective pair 
interaction between pseudo-macroions and an associated volume energy. The volume energy, which contributes to 
the total free energy, must be included when calculating thermodynamic properties of charged colloids modeled by 
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an effective pair potential |14|- |17| , |24| , p9[ . Noting the correspondences, microion «-> electron and macroion <-> metallic 
ion, the response approach is the colloidal equivalent of the widely-used pseudo-potential theory of metals (2^j||] . 

In previous work [ p0| , the response approach was extended to finite-sized macroions in deionized suspensions. This 
paper generalizes the theory to the case of arbitrary salt concentration, consistently taking into account (1) the 
volume excluded to the microions by the macroion hard cores and (2) the response of both counterions and salt ions 
to the macroion charges. The next section begins with a brief review of the response theory and then outlines our 



extensions of the theory. Section HI presents the main results - obtained within a linear response approximation 
- for an effective pair potential acting between pseudo-macroions and an associated volume energy, both of which 
consistently incorporate excluded volume effects. The influence of the volume energy on thermodynamic properties 
is illustrated by calculations of the osmotic pressure and bulk modulus. Comparisons with experimental data show 
that the counterions contribute a substantial fraction of the osmotic pressure and can soften or destabilize colloidal 



crystals. Finally, in Sec. IV we summarize and conclude. 



II. THEORY 
A. The Model 

Within the "primitive" model, the system comprises N m charged hard-sphere macroions of diameter a and charge 
— Ze (e being the elementary charge) and N c point counterions of charge ze suspended in an electrolyte solvent. Global 
charge neutrality constrains macroion and counterion numbers according to ZN m = zN c . Each macroion is assumed 
to carry a fixed charge, uniformly distributed over its surface. The solvent hosts N s pairs of salt ions in a uniform 
dielectric fluid characterized entirely by a dielectric constant e. For notational simplicity, we assume a symmetric 1:1 
electrolyte, consisting of N s point ions of charge ze and N s of charge —ze (i.e., same valence as counterions). The 
microions thus number N+ = N c + N s positive and iV_ = N s negative, for a total of = N c + 2N S . The system 
occupies a total volume V at temperature T and fixed salt concentration maintained by exchange of salt ions, through 
a semi-permeable membrane, with a salt reservoir. 

Denoting macroion and microion coordinates by {R} and {r}, respectively, the Hamiltonian of the system may be 
expressed in the general form 

H({K},{r}) = H m + H fl + H m+ + H in _. (1) 

The first two terms on the right side of Eq. (|l|) denote Hamiltonians for macroions and microions, respectively. 
Assuming the only relevant interactions to be steric and electrostatic, the bare macroion Hamiltonian, H m , is given 
by 



H m = K m + [vgs(\Ri - Rjl) + v mm (\TLi - R,Q] , (2) 



i,i=l 



A" m being the kinetic energy of the macroions, Uhs(|Rj — Rjl) the hard-sphere pair interaction between macroion 
cores, and v m m{r) — Z 2 e 2 /er the bare Coulomb interaction between a pair of macroions separated by center-to-center 
distance r > a . Similarly, the microion Hamiltonian takes the form 



Hp = 



*zZ v ns(\n - R/|) + - v ++ {\ri - rj\) 

1=1 j=l 



N+ N— _. N- 



E 



^^(Ir.-r.l) + I Yv-dri-rjl), (3) 



1=1 1 = 1 M=l 

where is the microion kinetic energy, ^Hs(l r i — Rj|) is the hard-sphere interaction between a point microion and 

a macroion core, and w++(r) = v (r) = — (r) = z 2 e 2 /er is the microion-microion Coulomb interaction. The last 

two terms in Eq. (|l|) are the macroion-microion electrostatic interaction energies, given by 

i=l 3=1 
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where v m ±(r) denotes the macroion-microion electrostatic pair interaction. For later reference, we note that Eq. 
also may be expressed in the form 



where 



H m± = J dr J dRp ± {r) Pm (R)v m ±(\r-R\), (5) 



N ± N m 
p±(r) = £«(r-r,), p m (R) = g 6(R - R,) (6) 

t=l j=l 

are the microion and macroion density operators, whose Fourier transforms are 

N ± N m 

P±(k) = ^exp(ik- r. ( ), p m (k) = ^ exp(ik • Ry). (7) 

Although v m ±(r) has the Coulomb form outside the macroion core radius, inside the core it has no unique definition. 
Thus, following van Roij and Hansen pj], we are free to choose v m ±(r) to be a constant for r < a/2 and take 

f^, r>a/2 

"- ±(r) M#«. —a 

where the parameter a will be specified (Sec. III.C) to ensure that the microion densities vanish within the core. 



B. Reduction to an Equivalent One-Component System 

With the Hamiltonian specified, we now turn to a statistical mechanical description of the system, our ultimate 
goal being the free energy. The canonical partition function is given by 

Z = ((cxp(-ff/fc B T)) Ai ) m , (9) 

the angular brackets symbolizing classical traces over microion and macroion degrees of freedom. Following standard 
treatments originating from the theory of simple metals |2^,^,^l|, we now reduce the two-component macroion- 
microion mixture to an equivalent one-component system by performing a restricted trace over microion coordinates, 
keeping macroion coordinates fixed. Thus, without approximation in this purely classical system, 

Z = (cxp(- H cS /k B T)) m , (10) 

where H c g = H m + F„ is the effective Hamiltonian of a one-component system of pseudo-macroions, and where 

F„ = - k B Thx{exp[-(H fl + H m+ + Hm-)/k B T])^ (11) 

may be physically interpreted as the free energy of a nonuniform gas of microions in the midst of macroions fixed at 
positions R^. Formally adding to and substracting from H the energy, Ey,, of a uniform background having a charge 



equal to that of the macroions, Eq. (11) may be recast in the form 



F M = -k B T\n( C xp[-(H'^ + H' m+ + H' m _)/k B T])^ (12) 

where H'^ = + and H' m ± — H m ± — E^/2. The advantage of this simple manipulation is that H'^ is the 
Hamiltonian of a classical, two-component plasma of microions, in a uniform compensating background, in the presence 
of neutral hard-sphere macroions. In order that the plasma be free of infinities associated with the long-range Coulomb 
interaction, the background must occupy the same volume as the microions. The background is thus excluded - along 
with the microions - from the macroion cores. The microions and background then jointly occupy a free volume 
V' = V(l — ij), which is the total volume reduced by the volume fraction of the macroion cores, 77 = (n /6)(N m /V)a 3 . 
The background energy is given explicitly by |31[ 
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1 o /" /" , z 2 e 2 /" (n + - n_)Zze 2 

£ b = -(n + -n_) 2 / dr / dr' - V / dr v 



2 V T y Jy, y v / e|r-r'| ^J v , e|r-R,| 
= -~(Ar + -AT_)(n + -n_)« ++ (0), (13) 

where n± = N±/V' = n± /(l — 77) are the effective mean densities of microions in the volume not occupied by the 

macroion cores and ?4 0) = N ± /V are the nominal mean densities. For later reference, we also define n s = n_ and 
n c = n + — rt_ as the effective densities of salt-ion pairs and counterions, respectively. In Eq. (fl3|), w_|_ + (0), defined by 

v ++ {0) = L dr — = (14) 

is the fc — * limit of the Fourier transform of v++(r). Although formally infinite, will be seen below to be 
identically cancelled by compensating infinities in and H m ± ■ 



C. Linear Response Approximation 

The theory presented thus far is exact, within the primitive model. The challenge remains to calculate the microion 
free energy [Eq. (|ll])]. One proposed strategy [H invokes density-functional theory to approximate F^, regarded as 
a functional of the microion densities, by performing a functional Taylor-series expansion about a uniform microion 
plasma. An alternative strategy (2^j2^], inspired by the pseudo-potential theory of metals, is to formally regard H' m ± 
as "external" potentials actin g up on a microion plasma and then approximate by perturbation theory. Following 
the second strategy, we write [glj 

F, = F p + [ dX{(H' m+ ) x + (H' m _) x ), (15) 
Jo 

where 

F p = -fc B rin(exp(-^;/fc B r)> /j (16) 

is the free energy of the reference microion plasma, occupying a volume V', in the presence of neutral hard-sphere 
macroions. The integral over A in Eq. ( |l5| ) corresponds physically to an adiabatic charging of the macroions from 
neutral to fully-charged spheres. The ensemble average ( ) A represents an average with respect to the distribution 
function of a system whose macroions carry a charge XZ . 

Further progress is facilitated by expressing (H{ n ±)\ in terms of Fourier components of the macroion and microion 
densities and of the macroion-microion interaction. From Eqs. (^|)-(Q), we have 

«±) A = ^X)^( fc )^±( k ))A^(- k ) + ^^^(*)<^( k ))A^(- k )]-^/ 2 ' 

Evidently {H' m ^. depends through /3±(k) upon the response of the microions to the macroion charge density. Re- 
garding the macroion charge as imposing an external potential on the microions, and assuming that the microion 
densities respond linearly to this potential, the Fourier components of the microion densities appearing in Eq. ( |T^ ) 
may be expressed in the form 

(p+(k)) A = \[ X++ (k)-X+-(k)}v m+ (k)p m (k), k^O, (18) 

and 

0Mk)) A = A [*+-(*) -X— (k)]v m+ (k)p m (k), fc^O, (19) 

where x±±(k) ar e the linear response functions of the reference two-component microion plasma and where we have 

used the symmetry relations Xh (k) — X h(^) an( i w ro +(fc) = — v m -(k). Note that for k = there is no response, 

since p±(0) — N±, as determined by the fixed numbers of microions. Substituting Eqs. ( ^8[ ) and ( |l9[ ) into Eq. ( |l7| ) 
and this in turn into Eq. (|l5|), and integrating over A, the microion free energy is given to second order in the 
macroion-microion interaction by 
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F» = F p + ^i 7 ^[ X++ (fc)-2 X+ _(fc) + X __(fc)][« m+ (fc)] 2 / 5 m (k) / 5 m (-k) 

+ N m (n + - n_) lim[v m+ (k)} - E h , (20) 
where again we have used the relation w m+ (fc) = —v m -(k). Correspondingly, the effective Hamiltonian takes the form 

1 Nm 1 

5 51 UHS (I R > ~ R iD + ^77 51 * mm ( fc ) [Pm(k)p m (-k) - Ay 



i N: 

H cS = K, 



2 ^ Vl J " 2V 

i,j=i k 



F p + 2^ E *(*) [WW] 2 /5m(k)/5 ro (-k) + AT ro (n + - n_) lim [WW] - ^b, (21) 

k^O 



where we have defined 



X {k) = x+ + (k)-2x+-(k)+x— (k). (22) 
Now rearranging terms, Eq. (|2l]) may be restructured and written in the formally simpler form 

1 Wm 1 
H eS = K m + 75 ^ "HsGRi - R-jl) + ttttt X/ " cff ( fc ) [/5m(k)/5 m (-k) - AT m ] + So 



2 ^ Jl/ 2F' 

«.J=i k 



- J] [VHs(|Ri-R,-|)+«efF(|Ri-Rj|)]+^, (23) 



2 



where 



v e s(r) = v mm {r) + v ind (r) (24) 



has the physical interpretation of an effective electrostatic pair potential between pseudo-macroions, which is the sum 
of the bare Coulomb potential and an induced potential whose Fourier transform is 

«tad(«0 = X(k)[v m+ (k)} 2 . (25) 

The final term in Eq. ( |23| ) is the volume energy, 

N r z i 

Eq = Fp + — Um u ind (r) + iV m (n+ - rz_) lim -— t>md(fc) + # m+ (fc) - £? b , (26) 

^ 7 — >0 fc^O L ZZj J 

which is a natural byproduct of the reduction to an equivalent one-component system. Although it has no explicit 
dependence on the macroion coordinates (see below), Eq evidently depends on the mean density of macroions and 
therefore can contribute significantly to the total free energy of the system. In passing, we note that the above 
expressions for the effective pair potential and the volume energy are analogous to expressions appearing in the 
pseudo-potential theory of metals |27],^,^5l],[32j if one substitutes for F p and x(fc), respectively, the energy and linear 
response function of the homogeneous electron gas (in the presence of a compensating background), and for v m +(k) 
the electron-ion pseudo-potential. 

Summarizing thus far, we have adopted the primitive model of charged colloids, formally reduced the macroion- 
microion mixture to an equivalent one-component system of pseudo-macroions, and applied a linear response approx- 
imation to the microion density, to obtain expressions for an effective electrostatic pair interaction [Eqs. ( p4| ) and 
(p5|)l and an associated volume energy [Eq. (pq)]. Practical calculations still require explicit specification of (1) the 
reference plasma free energy F p , (2) the plasma linear response functions x±±(^), and (3) the macroion-microion 
interaction v m +(k). In the next section we consider each of these in turn. 



III. RESULTS AND DISCUSSION 
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A. Reference Microion Plasma 



The free energy of the two-component reference plasma may be expressed as 

F p = F ic i + F COII + F cc — Eh, (27) 

where F^ and F COII are the ideal-gas and correlation contributions and F cc is the energy associated with Coulomb pair 
interactions between microions. It is important to emphasize that by associating the hard-sphere part of the total 
macroion-microion interaction with the microion Hamiltonian [Eq. (J3J) ] - required, since response theory does not 
apply to hard-sphere interactions - the reference microion plasma is implicitly restricted to the free volume outside 
of the macroion cores. As a consequence, the plasma is not strictly uniform, since the boundary conditions, imposed 
by the macroion surfaces, may induce nonuniformity. In general, the ideal-gas free energy is given by 



PEd = 



drpf (r) (ln[pf (r) A 3 ] - l) + J drpL 0) (r) (ln[pL 0) (r) A 3 ] - l) , (28) 



where (3 = l/fcrsT, A is the microion thermal de Broglie wavelength, and p± (r) are the nonuniform densities of 
positive and negative microions in an external field due to the macroion cores (but not their electric fields) . 

Now, for typical macroion charges and concentrations, counterion concentrations are in the /jM (10 -6 mol/1) range. 
If the salt concentration also falls in this range, then microion concentrations are low enough that the plasma is 
essentially uniform. In this case, F cc ~ E^ and the last two terms in Eq. ( p7| ) cancel each other. Furthermore, 
a plasma of such low concentration is weakly coupled, with coupling parameter T = z 2 e 2 / ek^Ta^ « 1, where 
a n = (3/47rn Al ) 1 / 3 is the microion sphere radius and = n + + n_ is the total microion number density. (In sharp 
contrast, electron plasmas in metals are typically characterized by T » 1.) The correlation free energy per microion 
then may be approximated by the Abe expansion |33j 

= _ _}_ rV2 + 0(r3)) (29) 

the leading term being the Debye-Huckel approximation [B4|. Thus, at low salt concentrations, if nonuniformities and 



correlations are ignored [[14 17 29], a reasonable approximation for the free energy of the microion plasma is 

f3F p ~ N+ [ln(n+A 3 ) - l] + N_ [ln(n_A 3 ) - l] 
Jfj, [ln(n M 

where x± — N±/N^ are the mean microion concentrations 



[ln(n p A 3 ) - 1 + ir + lniz: + + x_hixJ\ , (30) 



B. Linear Response Functions 

The linear response functions, Xij(k), hj = i, of the two-component reference plasma are simply proportional to 
the corresponding partial structure factors, Sij{k): 

Xij(k) = - fin^Sijik). (31) 

Liquid state theory [^5| now relates the partial structure factors to Fourier transforms of the pair correlation functions, 
hij(k), via 

Sij (k) — + XiXjn^hij(k). (32) 

The pair correlation functions are in turn related to Fourier transforms of the direct correlation functions, Cj,-(fc), by 
the Ornstein-Zernike (OZ) equation for mixtures: 

hij(k) = Cjj (k) + » M xicu (k)hij (fc), i,j,l = ± (33) 
l 

which serves in fact to define Cjj(fc). For such weakly-coupled plasmas as we encounter in charged colloids, the mean 
spherical approximation (MSA) provides a reasonable closure for the OZ equation. This amounts to approximating 
Cjj(r) by its asymptotic (r — > oo) limit, Cjj(r) ~ —f3vij(r), for all r, or equivalently, 



G 



dij (k) ~ - (ivy (k) 



4ir[3ziZje 2 

7k 2 : 



(34) 



where Zi, Zj — ±z. Since, in the MSA, c ++ (fc) = c (k) = — (fc) = c(k), it follows directly from Eq. (|33|) that 



h ++ (k) = &__(*) 



c(fe) 



1 — n^c{k) 



Substituting Eqs. (|32|), and (|35|) into Eq. (pi]) yields Xij(k), from which we obtain 



X++(k) - X+-(k) 



/3n H 



1 - n^c(k) 



j3n + 



l+K 2 /k 2 ' 



(35) 



(36) 



and 



where 



x+-(fc) - X — (fc) 



x(*) 



/3n_ 



/3n_ 



l-n M c(fc) l + K 2 /fc 2 ' 



/3n, 



1 - n fi c(A;) 



l + K 2 /fc 2: 



/4 7 rn AI z 2 e 2 \ 1/2 / 47r4 0) z 2 e 2 



1/2 



(1 - 7?)£A:bT 



(37) 



(38) 



(39) 



and 



,(0) 



Np/V = 7i M (l — rf) is the total nominal microion number density. As will be seen below, the parameter k 



plays the role of the Debye screening constant (inverse screening length) in the microion density profiles and in the 
effective pair interaction. 



C. Microion Density Profiles 



Specifying the macroion-microion interaction amounts to determining the value of the parameter a in Eq. (pfo that 
ensures vanishing microion densities inside the macroion cores. This in turn requires a calculation of the real-space 
microion density profiles. The first step of this calculation is to Fourier transform Eq. (||), with the result 



Vm±(k) 



inZze 2 
ek 2 



(1 — a) cos(fca/2) + a 



sm(fcer/2) 
ka/2 



(40) 



Now substituting Eqs. (36), (|37|), and ( [4Gj ) into Eqs. (18) and Eqs. ( |19| ) gives for the k ^ Fourier components of the 
microion densities 



p±(k) = ± x± — 



Z \k 2 + K 2 



(1 — a) cos(fccr/2) + a 



sin(fccr/2) 
ka/2 



N n 



^ exp(zk • Rj), k ^ 0, 



(41) 



where the sum is over the positions Rj of the macroions. Next inverse transforming Eq. (^), while respecting the 
k — > limits, p±(Q) — N±, we obtain 



P±(r) 



± z ± £>>(|r-R,-|), \r-Rj\>a/2 



(42) 



ar±$><(k-:Ril) 5 |r-Ri|<a/2, 



where = x_n + + x + n^ is the bulk density of positive or negative microions (far from any macroion). Note that 
in general p^ ^ n s , although p^ — > n s in the limit x± — > 1/2 (n c /n s — > 0). In Eq. (|4^), p>(r) and p < (r) are 
single-macroion orbitals, given by 
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P>(r) 



Z k 2 \, s , , s sinh(/tcr/2) 
(1-a) cosh(«o-/2) + a 



z 4tt 



na/2 



exp(— nr) 



>a/2, 



and 



z An 



-1 + a- 



exp(-Kcr/2) — '-, 



na/2 



r < a/2. 



Vanishing of p<(r) for r < a/2 is evidently ensured by setting 

na/2 



o 



1 + na/2' 

Finally substituting this choice for a back into Eq. (E3J) specifies the r > a/2 orbital as 



Z k 2 exp(na/2) exp(— nr) 
P> ^> = ~An~ 1 + na/2 r ' 



r > a/2, 



(43) 



(44) 



(45) 



(46) 



which is automatically normalized to the correct number of counterions per macroion (Z/ z). The corresponding 
microion density profiles are the linear combinations 



. Z k exp(Kcr/2) exp(— «|r — Rd) 

P±[T) = Poo ± ar± — 7- x 



z An 1 + kct/2 



E 



|r-Rj| >a/2. 



(47) 



Expression (fig ) is seen to be of precisely the same form as the Debye-Hiickel expression for the density of electrolyte 
ions around a macroion [Q. A significant distinction lies, however, in the definition of the screening constant, n. 
Whereas the Debye-Hiickel k depends on the nominal bulk density of electrolyte ions, our k [Eq. ([59])] depends rather 
on the effective mean microion density (in the volume unoccupied by macroions). The importance of redefining the 
usual k in this way, particularly for concentrated suspensions, has been noted previously by Russel and coworkers [p6[. 



D. Effective Pair Interaction and Volume Energy 



We are now in a position to derive the main results of the paper. Considering first the effective electrostatic 
pair interaction between pseudo-macroions, we proceed by substituting Eq. into Eq. (|t)]), obtaining for the 
macroion-microion interaction 



i> m ±(k) 



AnZze 2 



ek 2 



1 



1 + na/2 



cos(ka/2) + k 



sin(ka/2) 



Next substituting Eqs. ( p8|) and (|48| ) into Eq. (25) yields the induced potential: 



Wind(fc) 



2nZ 2 e 2 



1 



1 + cos(fcer) + 2k 



sin(fcer) 



, f — cos(fcer) 



ek 2 \l + na/2) \k 2 + 

Fourier transformation of Eq. (^9|) is a straightforward calculation, with the result 

2 

exp 



(48) 



(49) 



Z 2 e 2 ( exp(Kcr/2) \ exp(-Kr) _ Z 2 e 2 
£ \ 1 + na/2 ) r er > 

Z^e 



(1+^72) [( 2 + ^)--5« 2 



r > a 



r < a. 



(50) 



Finally substituting Eq. (M) into Eq. (EJ) , we obtain an explicit expression for the effective electrostatic pair potential: 



v e s(r) 



Z 2 e 2 ( exp(Kcr/2)^ exp(— kt) 
e V 1 + K(7 / 2 



r > a. 



(51) 



This result is seen to be identical in form to the electrostatic part of the DLVO effective pair potential || , which is 
usually derived by linearizing the Poisson-Boltzmann equation. The only distinction between our pair potential and 
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the DLVO potential lies in the definition of the screening constant, ours [Eq. (|3^)] being a factor (1 — r\)~ x l 2 larger 
than the usual DLVO k to account for exclusion of microions from the macroion cores. 

Now the volume energy may be explicitly determined from Eq. (^6|). It follows immediately from Eq. ([SO]) that 

Z 2 e 2 k . . 

hmv ind {r) = — , 52 

r^O e 1 + KG j 1 



from Eq. (49) that 



- ,,s f Z Y~ *Z 2 e 2 a 2 l + Ka/6 4irZ 2 e 2 

Lmwind fc = - - «++(0) + — + — , 53 

fe^o \ z J e 1 + kcx/2 en z 



limi, m+ fc - --«++(0) + — ^. 54 

fe^o z 2e 1 + Kcr/2 



and from Eq. (pf) that 



Substituting Eqs. ( p2J ) - (|5J) into Eq. (f26|), and using approximation fl30|), we obtain for the volume energy 

PE = N+ ln(n+A 3 ) + iV_ ln(n_A 3 ) - N m ^l K -— - I C^Z^l! (55) 

' + 1 + 7 1 7 2e I + kct/2 2 7V+ + iV_ ' 17 

neglecting irrelevant constants. Note that the infinities associated with the k — > limits formally cancel one another, 
as t hey m ust p7|. The first term on the right side of Eq. (|55|) is the ideal-gas plasma free energy, discussed in 



Sec. HI A. The second term, which accounts for the electrostatic energy of interaction between the macroions and 
their screening clouds of counterions, is equivalent to one half the interaction energy were all the counterions to be 
placed at a radial distance n^ 1 from the centers of their respective macroions. The final term corresponds to the 



k — > limit in Eq. (|26|). Our result for the volume energy is very similar to that derived by van Roij et al. []14 15 
from a density-functional expansion, differing only in the manner in which exclusion of microions from the macroion 
cores is incorporated. While Eq. (|5^) incorporates excluded volume effects through a dependence of the screening 
constant [Eq. (|39|)1 on the effective microion density, van Roij et al. incorporate them through an additional term in 
the volume energy [Eq. (61) in Ref. [p^[]. 

In closing this section, we remark on the range of validity of the theory. First, although the linear response 
approximation presupposes relatively weak microion response to the macroions, and thus weak screening, the general 
form of the screened- Coulomb pair potential in bulk suspensions is broadly supported by Poisson-Boltzmann cell 
model calculations |To|| , ab initio simulations Jl^,^3|, and experiments Q. Second, the excluded volume corrections 
incorporated in the modified screening constant, k, may become significant even in the weak-screening regime for 
concentrated suspensions of weakly charged macroions. Finally, although the theory neglects (in mean-field fashion) 
fluctuations and correlations in the microion densities, Monte Carlo simulations and cell model calculations [E3 for 
spherical macroions suggest that such correlations contribute only marginally to the total free energy. 



E. Osmotic Pressure and Bulk Modulus 



Being independent of the macroion coordinates, the microion volume energy, Eq, appears simply as an additive 
term in the total Helmholtz free energy of the system: F = F m + Eq, where F m is the free energy of the equivalent 
one-component system of pseudo- macroions interacting via the effective pair potential v e s(r). Correspondingly, any 
thermodynamic quantity derived from this free energy may be decomposed into effective macroion and microion 
contributions. Since Eq [Eq. (|5q)1 depends on the mean macroion density - both explicitly and implicitly through k - 
it can significantly influence thermodynamic properties of the system, especially at low salt concentrations [p5|- p7| , p0[ . 

As illustrations, we consider the osmotic pressure and bulk modulus. A colloidal suspension in equilibrium, through 
a semi-permeable membrane, with a reservoir of salt solution exerts an osmotic pressure, II = P — P r , defined as the 
difference between the pressure of the system, P, and that of the reservoir, P r . Treating the reservoir as an ideal gas 
of N r salt ion pairs in a volume V T , we have /3P r = 2N T /V T . Chemical equilibrium is characterized by equality of the 
chemical potentials of salt ion species exchanged between the system and the reservoir. The chemical potential of the 
salt, defined as the change in free energy upon adding a salt ion, includes a contribution arising from the effect of 
salt concentration on the macroion-macroion interaction (through k). Thus, in general, the salt concentrations of the 
system and reservoir are nontrivially related. However, for systems sufficiently dilute that the macroion contribution 
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may be ignored - an assumption we make here - the condition for chemical equilibrium may be approximated by 
equality of the reservoir salt density, N x /V x , and the effective salt density of the system, n s — (N s /V)/(1 — 77). Note 
that the effective salt density exceeds the nominal salt density, N s /V, by the ratio of the total volume to the free 
volume unoccupied by the macroion cores. The distinction here between nominal and reservoir salt densities is akin 



to that between nominal and reservoir polymer densities in colloid-polymer mixtures |38|. The reservoir pressure is 
then given by 

I3P X ~ 2n s . (56) 
The total pressure (or equation of state) of the system, P = P m + Pq, comprises a macroion contribution, P m , and a 



microion contribution, Pq = —(dEo/dV)N m ,N a - Combining Eqs. (pq) and (56), we obtain 



167T z A (1 + n<7/2y 



The same result is obtained for arbitrary macroion concentration in the limit of zero salt concentration (n s — > 0), 
in which case P x = 0. The second and third terms in Eq. fl57| ) represent, respectively, the ideal-gas pressure of the 
counterfoils and a van der Waals-like adjustment that accounts for the attraction between counterions and macroions. 
With increasing counterion density, these two terms compete with each other, the attractive term acting to reduce 
the total osmotic pressure. For weak microion screening (na < 1), where the electrostatic fields of the macroions 
are relatively weak and the microion densities close to uniform, Eq. ( pT|) should reasonably approximate the osmotic 
pressure. In fact, in the limit no — > 0, our result naturally tends to the correct ideal-gas limit, (5Pq — > n c . For 
stronger screening (kg > 1), where the microion densities are more nonuniform, nonlinear response effects may 
become significant. 

The bulk modulus (or inverse compressibility), defined by B = — V(dIL/dV) N m ,N a , may be similarly expressed in 
the form B = B m + Bq, where B m and Bq are the macroion and microion contributions, respectively. From Eq. (|57]), 
we immediately obtain 

o Rn 3 _ o B * , "c<x 3 3 Z 1 (^) 3 (1 + W6) , . 



which again includes repulsive and attractive microion terms. In Eqs. (p7j) and (p8|), the macroion contributions, 
P m and P m , are understood to be obtained from a theory (or simulation) of a one-component system of particles 
interacting via the effective pair potential [Eqs. ( p9[ ) and (pi])]. In practice, the macroion charge, notoriously difficult 
to extract from experiment, is usually replaced by an adjustable parameter, the effective or renormalized charge, 

Z* . 

As a test of our results, we compare, in Fig. 1, the osmotic pressure predicted from Eq. ( |57| ) with the recent 
experimental measurements of Reus et al. [[59) for a colloidal fee crystal in a highly deionized (n s = 0) aqueous solvent 
at room temperature (Bjerrum length Ab = /3z 2 e 2 /e = 0.714 nm). As an approximation for the macroion pressure, 
P m , we use results of integral-equation calculations based on the virial equation with HNC closure for the liquid-state 
pair distribution function For the effective macroion charge, we take the value Z* = 700 estimated by Reus 

et al. to best match their phase diagram to the simulations of Robbins et al. [ fi"H . 

The microion contribution is seen to make the dominant contribution to the total osmotic pressure and to sub- 
stantially improve the agreement between the one-component model and experiment, particularly at lower volume 
fractions (77 < 0.07). The last term in Eq. ([57]) clearly is essential to reduce the rapidly increasing counterion ideal-gas 
pressure. The discrepancies at higher volume fractions (r/ > 0.07) might be attributed, at least partially, to an under- 
estimate of P m by the liquid-state theory. They may also reflect nonlinear response of the counterions and associated 
effective many-body interactions between pseudo-macroions. Future work will address influences of effective triplet 
interactions on the osmotic pressure 0]. It should be mentioned that the Poisson-Boltzmann cell model [^2|39| 
(upper curve in Fig. 1) matches the experimental data well, especially at higher rj. However, while cell models, which 
consider the distribution of microions within a Wigner-Seitz cell centered on a single macroion, are limited to periodic 
crystals, the more general one-component model applies to any thermodynamic phase. 

A more stringent test of the theory is presented by the bulk modulus - the curvature, with respect to density, of 
the free energy density. In recent experiments, Weiss et al. [[l3| determined the bulk modulus of colloidal fee crystals 
suspended in a deionized, aqueous solvent at room temperature (Ab = 0.714 nm) by measuring the long-wavelength 
limit of the static structure factor. For two samples, distinguished by nearest-neighbor distances, a = (3/47rn m ) 1 / 3 = 
2.5 /im and 3.25 /mi, the measured bulk moduli were argued to be lower than the predictions of DLVO theory, as 
estimated on the basis of an approximate elastic theory for the macroion contribution, P m . For the denser crystal, 
the measured value was B = 0.016 ± 0.005 Pa, less than a third of the estimated "DLVO" value of B = 0.052 ± 0.005 
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Pa. This analysis ignores, however, the counterion contribution associated with the volume energy. Figures 2 and 3 
present predictions, computed from Eq. (|58|), for the counterion contribution, Bq. These results demonstrate that for 
sufficiently high effective macroion charge and volume fraction the counterion contribution may become negative. It is 
essential to include this contribution in the total bulk modulus before comparing the DLVO theory with experiment. 
In Fig. 2, the cross-over point at Z* ~ 7100 may be compared with the effective charges, Z* ~ 6100 for a = 2.5 
/im and Z* ~ 5200 for a — 3.25 fiiii, estimated by Weiss et al. for isolated pairs of spheres in the infinite dilution 
limit. However, lacking reliable knowledge of Z* in the crystal phase, we forgo here a more quantitative analysis. 
The qualitative message is nevertheless clear: at sufficient concentration, the counterions may act to lower the bulk 
modulus, softening or even destabilizing (B < 0) the crystal. 



IV. CONCLUSIONS 



In summary, by reducing a model colloidal suspension of charged hard-sphere macroions and point microions to 
an equivalent one-component system, and approximating the microion response to the macroion charge using linear 
response (second-order perturbation) theory, we have derived an effective electrostatic pair interaction [Eq. (|5l])] and 
an associated microion volume energy [Eq. fl5J)]. The volume energy, which depends on the average macroion density, 
accounts for both the microion entropy and the macroion-microion interaction energy. The effective interaction, which 
governs the dynamics of the macroions, is of precisely the conventional DLVO form for finite-sized macroions, but 
incorporates excluded volume corrections through the dependence of the screening constant on the effective density 
of microions in the free volume between macroion cores. 

The total free energy of the system is the sum of the volume energy and the free energy of the equivalent one- 
component system of pseudo-macroions. From the free energy, we have derived simple analytic expressions for the 
osmotic pressure and bulk modulus. Comparison of theoretical predictions with experimental data for deionized 
suspensions of highly charged macroions shows that the microions can significantly contribute to the thermodynamic 
properties, beyond their role in screening the bare Coulomb interaction between macroions. In particular, the volume 
energy largely accounts for the observed magnitude of the osmotic pressure and qualitatively explains measurements 
of bulk modulus lower than predicted by the conventional one-component DLVO theory. Several recent studies have 
predicted similar influences of volume energies on the phase behavior of charged colloids p^- jl7| , p^ | . 

The theory presented here can be straightforwardly generalized to include nonlinear response of microions f42]] and 
thereby used to assess the relative importance of effective many-body interactions |l||23|,^5| and associated corrections 
to the effective pair potential and the volume energy. Related applications are to colloid-surface interactions and to 
interactions between colloids in the vicinity of a surface, which experiment Q and theory Jl8[ suggest may become 
attractive. Work along these lines is in progress. 
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Figure 1 




FIG. 1. Osmotic pressure II vs. macroion volume fraction r) for an fee crystal of spherical macroions (diameter a = 102 
nm) suspended in a salt-free, aqueous solvent at room temperature (Ab = 0.714 nm). Symbols: experimental data of Reus et 
al. |}9|; solid curve: prediction of linear response theory [Eq. (p7[)] with effective macroion charge Z* = 700 and HNC-virial 
pressure for macroions (see text); long-dashed curve: counterion contribution IIo; short-dashed curve: Poisson-Boltzmann cell 
model prediction p2[p9|. Over this range of volume fractions, the screening constant increases from zero to kg ~ 4. 
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FIG. 2. Counterion contribution to bulk modulus Bo vs. effective macroion charge Z* for fee crystals of spherical macroions 
(diameter a — 654 nm) suspended in a salt-free, aqueous solvent at room temperature. Parameters are chosen for comparison 
with ref. Solid curve: nearest-neighbor distance a — 2.5 fim (macroion volume fraction r\ = 0.0133); dashed curve: 

a = 3.25 /im (rj = 0.00475). For Z* > 7100 the counterion contribution is negative. The arrow indicates the estimated 
macroion contribution to the bulk modulus of the denser crystal (see text). The maximum screening constant is kcj — 1.9. 




FIG. 3. Counterion contribution to bulk modulus Bo vs. macroion volume fraction r\ for same system as in Fig. 2, but 
for three different effective macroion charges. Solid curve: Z* = 6000; long-dashed curve: Z* — 7000; short-dashed curve: 
Z* — 8000. For sufficiently high Z* , the counterion contribution may be negative over a significant range of volume fractions. 
The arrow indicates the estimated macroion contribution to the bulk modulus of the denser crystal jl3|. The maximum 
screening constant is kct = 2.5. 
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